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The conformal covariance of correlation functions is checked in the second-order transition in- 
duced by random bonds in the two-dimensional eight-state Potts model. The decay of correlations 
is obtained via transfer matrix calculations in a cylinder geometry, and large-scale Monte Carlo 
simulations provide access to the correlations and the profiles inside a square with free or fixed 
boundary conditions. In both geometries, conformal transformations constrain the form of the spa- 
tial dependence, leading to accurate determinations of the order parameter scaling index, in good 
agreement with previous independent determinations obtained through standard techniques. The 
energy density exponent is also computed. 
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It is well known that quenched randomness can deeply 
affect the critical properties at second-order phase tran- 
sitions [1], and is liable to smooth first-order transitions, 
eventually leading to continuous transitions [2,3]. 

In random systems, owing to strong inhomogeneities 
inherent in disorder, the usual symmetry properties re- 
quired by conformal invariance (CI) [4] do not hold. How- 
ever, by averaging over disorder realizations (denoted 
by [. • -]av)i translation and rotation invariance are re- 
stored, and it is generally believed that conformal invari- 
ance techniques should apply in principle. Recent results 
based on this assumption have recently been obtained at 
randomness-induced second-order transitions [5,6], but 
clear evidence of the validity of conformal invariance is 
still missing. The question is of both fundamental and 
practical interest. From the fundamental point of view, 
situations are known where a diverging correlation length 
does not guarantee the validity of CI. A few years ago, 
lattice animals were indeed found to be not conformally 
invariant although they display isotropic critical behav- 
ior with correlation lengths satisfying the usual scaling 
£ ~ L with the system size [7]. If conformal invariance 
works, on the other hand, its powerful techniques might 
be applied with no restriction to investigate the critical 
behavior of 2D random systems [8]. 

In the 2D random-bond Ising Model, both analytic [9] 
and numerical results [10] are available. Transfer ma- 
trix (TM) calculations were also used to study the cor- 
relation function decay along strips and to compute the 
conformal anomaly (defined as a universal amplitude in 
finite-size corrections to the free energy) [11] and, since 
disorder is marginally irrelevant in the 2D Ising model, 
conformal invariance techniques were indeed efficient. At 
randomness- induced second-order phase transitions, a di- 
rect comparison between the results deduced from confor- 
mal invariance, and standard techniques, such as finite- 
size scaling (FSS) Monte Carlo (MC) simulations, have 
nevertheless not yet been made. After the pioneering 
work of Imry and Wortis, the first large-scale MC simu- 
lations devoted to the influence of quenched randomness, 



in a system whose pure version undergoes a strong first- 
order phase transition, were applied to the case of the 
eight-state random-bond Potts model (RBPM) [12]. 

The Hamiltonian of the system with quenched ran- 
dom nearest-neighbor interactions is written: —pTl = 
S(rr') Krr'fia r ,a r j ) where the spins, located at sites r of 
a square lattice, take the values cr r = 1, 2, . . . , q. The 
coupling strengths are allowed to take two different val- 
ues K and K' = Kr with probabilities p and 1 — p, re- 
spectively. If both couplings are distributed with the 
same probability, p = 0.5, the system is, on average, 
self-dual and the critical point is exactly given by dual- 
ity relations: {e K - - l){e K " r - 1) = q [13]. This model 
has again been carefully investigated recently using both 
MC simulations [14,15] and TM calculations [16]. The 
bulk magnetization scaling index x a = (3/v is q de- 
pendent [16]. Up to now the most refined direct esti- 
mate of x h a in the case q — 8 is probably a FSS anal- 
ysis due to Picco {x h a — 0.150 — 0.155) [15], where it 
was shown that the random fixed point is reached in 
the range r = 8 — 20, while crossover effects perturb 
the results outside this domain. As mentioned above, 
TM calculations have already been performed on the 
RBPM [16], but with a rather weak disorder, r = 2 lead- 
ing to x b a — 0.142(1), while standard FSS results are in 
the range x b a = 0.158 — 0.175 [15]. The discrepancy may 
presumably be attributed to crossover effects. 

Due to the increasing literature devoted to second- 
order induced phase transitions, in this Rapid Commu- 
nication, our investigation deals with the self-dual eight- 
state RBPM. We report results of TM calculations on the 
strip and MC simulations in a square geometry, which 
support the assumption of conformal covariance of or- 
der parameter correlation functions and profiles. Wc 
particularly compare independent determinations of the 
bulk magnetization exponent, resulting in good agree- 
ment between standard and conformal invariance tech- 
niques. Different large-scale simulations are performed 
in two types of restricted geometries (strips and squares) 
with a ratio r chosen in the range r = 8 — 20. The 
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crossover regime r = 2 is also investigated. 

In the strip geometry, we used the connectivity trans- 
fer matrix formalism of Blote and Nightingale [17]. For 
a strip of size L with periodic boundary conditions, the 
leading Lyapunov exponent follows from the Furstenberg 
method: A (L) = ]im m ^ ao ^ In || (IlfeLi T fc) l u o)||, where 
Tfc is the transfer operator between columns k — 1 and 
k, and | Vq) is a suitable unit initial vector. The leading 
Lyapunov exponent determines the free energy: fo{L) = 
— L~ 1 Aq(L). For a specific disorder realization, the spin- 



spin correlation function (G a (u)} 



, where 



(. . .) denotes the thermal average, is given by the prob- 
ability that the spins along some row, at columns j and 
j + u, are in the same state: 



<oig,(nitr lT, fe )d J+M |o) 



(oinitr 



'T fc |0> 



(1) 



where 1 0) is the ground state eigenvector, T' k is the trans- 
fer matrix in the extended Hilbert space [18], gj is an op- 
erator that identifies the cluster containing o~j, and dj+ u 
gives the appropriate weight depending on whether or 
not (Jj+u is in the same state as o~j . 

It is well known that in disordered spin systems, the 
strong fluctuations from sample to sample can induce av- 
erage difficulties [19]. For that reason we paid attention 
to check, by a careful analysis of the correlation function 
probability distribution, that self-averaging problems do 
not alter the mean values [20] . In order to reduce sample 
fluctuations, we furthermore considered canonical disor- 
der, a situation in which exactly the same amount of 
both couplings is distributed over the bonds of the sys- 
tem. The computations are then performed with 10 6 
iterations of the transfer matrix, and the final Lyapunov 
exponent is averaged over 20 disorder configurations. 



Since our purpose is to check the predictions of confor- 
mal symmetry, which are supposed to be satisfied by av- 
erage quantities, i.e., [(G a (u))] av , our first aim is to show 
that, in spite of the lack of self-averaging, our numerical 
experiments lead to well-defined averages. In Ref. [16], 
Jacobsen and Cardy argued that in the RBPM, InG is 
self-averaging while G is not. Exploiting duality, they 
computed [ln(G a (u))} a v via the free energy of a system 
in the presence of a seam of frustrated bonds and a cu- 
mulant expansion enabled them to deduce the behavior 
of [(G a (u))} av . In our approach, we are first interested in 
the probability distribution of the correlation function. 
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FIG. 1. Average correlation function, most probable value, 
and fourth order cumulant expansion obtained from 63 436 
replicas for a strip of size L — 6 (r = 10). 

The most probable value G™ p (u) and the average cor- 
relation function [{G a {u))] a vi as well as the averaged log- 
arithm [ln(G cr (u))] au , can then be deduced at any value 
of u. Compatible behaviors are found for G™ p (u) and 
e [ ln (G„ (u))]av ; which confirm the essentially log-normal 
character of the probability distribution [20], in agree- 
ment with Cardy and Jacobsen. It is thus necessary to 
perform averages over a larger number of samples for 
[(Gv(u))]av than for [ln(G cr (u))] al , to get the same rel- 
ative errors. A cumulant expansion enables us to re- 
construct [(G a (u))] av and to compare to the values ob- 
tained by averaging directly over the samples. The re- 
sults in Fig. 1 strengthen the credibility of the direct 
average and also clearly show that the cumulant expan- 
sion up to fourth order still strongly fluctuates compared 
to [(G a (u))] av . 

We will now concentrate on the results that follow 
from the assumption of conformal covariance of the aver- 
aged correlation functions. In the infinite complex plane 
z = x + iy at the critical point, the correlation func- 
tion exhibits the usual algebraic decay [{G a {r))] a v — 
const xr" 2l », where r =\z\. Under a conformal mapping 
w(z), the correlation functions of a conformally invariant 
2D-system transforms according to 

G a {w 1 ,w 2 )^\w'{z 1 )\-<\w'{z 2 )\-< G a ( Zl ,z 2 ). (2) 

The logarithmic tranformation w — -k- In z is known to 
map the z plane onto an infinite strip w = u + iv of 
width L with periodic boundary conditions in the trans- 
verse direction. Applying Eq. (2) in the random sys- 
tem, one gets the usual exponential decay along the strip 
[(G CT (u))] a „ = const x exp (— =£-x b a v), where the scaling 
index x\ can be deduced from a semilog plot. For each 
strip size (L — 2 — 9), we realized 40 000 disorder con- 
figurations in four independent runs, which allowed us to 
define mean values and error bars. The exponent follows 

from an exponential fit in the range u = 5 10L. For 

r = 10, the resulting values plotted against L~ x , con- 
verge in the L — > oo limit 
x b a = 0.1496 ±0.0009 shown in Fig. 2 



is obtained in the case r 



towards the final estimate 
A consistent value 
20, while in the weak disor- 
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der limit the behavior is drastically different, and strong 
crossover effects would be expected for larger strip sizes. 



0.18 



0.16 



0.14 



0.12 



0.0 



r=10 [0.1496(9)] 
r=20 [0.1474(19)] 
r=2 [0.1256(5)] 




0.1 



0.2 0.3 

1/L 



0.4 



0.5 



10 



A 

V 10"' 




0.155(7) 

0.153(8) 

0.152(9) 

0.151(8) 

0.150(10) 

0.150(8) 



10 • 



10 

K(Q 



10" 



FIG. 3. Rescaled correlation function along six u = const 
curves in the square (r = 10). 



FIG. 2. Magnetic scaling index deduced from the algebraic 
decay of the average correlation function along the strip of 
size L as a function of L" 1 . 

Another restricted geometry, already used in pre- 
vious FSS MC simulations, is the square geometry. 
The Schwarz-Christoffel mapping ( — 4F(z, k), z = 



sn 



=«n f»££ 



sn ( > where F(z,k) is the elliptic integral 

of the first kind and sn(£, k) the Jacobian elliptic sine [21] , 
transforms the upper half-plane inside the interior of a 
square -a/2 < Re(C) < a/2, < Im(C) < a. Here, 
K = K(k) is the complete elliptic integral of the first kind 
and the modulus k is a solution of K(fc)/K(Vl — k 2 ) = \. 
The correlation function in the semi-infinite geometry is 
known to take the form [22] 



G a {z 1 ,z) = Cst.{y 1 y)- x ^{u), 



(3) 



where the dependence on u> 



VlV 



of the universal 



scaling function ip is constrained by the special confor- 
mal transformation. In the random situation one can 
again use the transformation equation (2) to write the 
correlations between £i = i, close to a side of the square, 
and any point inside it, as follows: 

[{Ga(0)U = const x [| C(z) | Im(z(C))r^M- (4) 

Taking the logarithm of both sides, the bulk critical ex- 
ponent x b a can thus be deduced from a linear fit along 
u> = const curves in the square: 

ln[{G a (0)]av =const / -x*lnK(C)+lnV'(w), (5) 

with k(C) = Im(z(C))| [1 - z 2 (()] [1 - fcV(C)] |" 1/2 . 

The results are shown in Fig. 3 (r = 10), in which the 
Swendsen-Wang algorithm has been used [23] for systems 
of size 101 x 101 and an average was performed over 3000 
disorder realizations. Averaging the results at different 
ui's, one obtains x h a — 0.152 ±0.003. Again, a compatible 
value is found at r = 20 while it disagrees at r = 2. 



Owing to the unknown scaling function the de- 

termination is not extremely accurate, since a few points 
are used for the fits. It can nevertheless be improved if 
one considers the magnetization profile inside a square 
with fixed boundary conditions. Since it is a one-point 
function, its decay from the distance to the surface in 
the semi-infinite geometry is fixed, up to a constant pref- 
actor {(a(z))] av ~ y~ x " ■ The local order parameter is 
defined, according to Ref. [24], as the probability for the 
spin at site £ in the square, to belong to the majority 
orientation. After the Schwarz-Christoffel mapping, one 
gets the following expression for the average profile in the 
square geometry: 



const x 



Im(z(C)) 



(G) 



This expression, of the form [(c(C))]an = [f( z )/y] x ", 
holds for any point inside the square. It allows an 
accurate determination of the critical exponent {x a = 
0.1499 ± 0.0001 for r = 10) via a log-log plot shown in 
Fig. 4. We note that in the case r = 2, the correspond- 
ing curve exhibits a crossover between small and large 
distances, with clearly different exponents close to 0.128 
and 0.160, respectively. 
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FIG. 4. Rescaled magnetization and energy density profiles 
inside the square for 3000 disorder realizations (r = 10). The 
power law fits are over 100 2 data points. 
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TABLE I. Comparison between FSS and CI determina- 
tions of the bulk magnetic scaling dimension in the q — 8 
RBPM. The quantity that was studied is indicated in the ta- 
ble as well as the geometry and the numerical technique. 



FSS(MC) 




Conf. Inv. 


square 


strip 


strip square 


SW a W b 


TM C 


TM d SW 


[{M b )\ [{M b )\ 


[<G(u)>] 


l(G(u))] [(G(CW [M0)] c 



r = 10 0.153(1) 

0.153(3) 0.152(2) 0.1496(9) 0.152(3) 0.1499(1) 
r = 20 0.150(1) 

0.152(2) 0.1474(19) 0.146(7) 0.1462(2) 

r = 2 0.167(2) 0.142(1) 

0.158(3) 0.142(4) 0.1256(5) 0.283(21) 0.128-0.160 

a MC simulations (Swendsen-Wang algorithm, a < 100, ~ 500 
samples) from Ref. [14]. 

b MC simulations (Wolff algorithm, ~ 10 5 samples) from Ref. 
[15]. The two values refer to fits in the ranges a = 20 — 100 
and 50-200. 

C TM calculations (L = 1 — 7, 10 2 samples) from Refs. [6,16]. 
d TM calculations (L = 2 - 9, 4 x 10 4 samples), this work. 
e MC simulations (a = 101, 3 x 10 3 samples), this work. 

A summary of our results, compared to independent 
FSS determinations of the magnetic scaling index, are 
given in Table I. Apart from the crossover regime at 
r = 2, the agreement is quite good and clearly in favor 
of the validity of the assumption of conformal covari- 
ance of correlation functions and profiles at the random 
fixed point. On the other hand, in the crossover regime, 
in which the fixed point is not yet reached, the different 
techniques lead to different results, and none of them has 
a real meaning, since they would presumably be affected 
by strong crossover effects at very large sizes. Since the 
study of the critical profile with fixed boundary condi- 
tions inside the square is very accurate, the energy scal- 
ing index x\ = d — 1/u can also be obtained through the 
same technique. This dimension has been calculated in 
previous studies [6,14] but different results were obtained, 
possibly contradicting the inequality x\ > 1. Here, the 
local energy density [(s(())]av is computed by the av- 
erage over four bonds on each plaquette. It includes a 
constant bulk contribution [(eo)]at>, which is obtained by 
the extrapolation to y — > oo of the profiles with free and 
fixed boundary conditions (BC), respectively. At r = 10, 
both BC's yield consistent values [(eo)] at > = 0.6974 and 
0.6978. The quantity A[(e(())] av = [(s(())] av - 0.6976 is 
then studied as in Eq. (6), and leads to x h e = 1.011±0.008 
(Fig. 4). Since the numerical data are very small, there 
is an important dispersion, but the accuracy of the fit 
remains correct due to the huge number of data points 
(100 2 ). 

In this Rapid Communication, we have shown that 
conformal invariance techniques can be successfully ap- 
plied to random systems, and that they lead to refined 
investigations of the critical properties. The accuracy, 
compared to standard techniques, is increased, especially 



through the magnetization profile inside a square where 
all of the lattice points enter the fit. The critical ex- 
ponents are finally quite close to the rational values 
x\ = 3/20 and x\ = 1. 
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